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Abstract 

n 

= time step index, shedding frequency 

Unsteady rotor wake interactions with the empennage, tail 
boom, and other aerodynamic surfaces have a significant 

P 

= transformation parameter 

influence on the aerodynamic performance of the 
helicopter, ride quality, and vibration. A Computational 

P 

= pressure 

Fluid Dynamic (CFD) method for computing the 
aerodynamic interaction between an interacting vortex 

Q 

= transformation parameter 

wake and the viscous flow about arbitrary two- 
dimensional bodies has been developed to address this 

r 

= radial variable 

helicopter problem. The vorticity and flow field velocities 
are calculated on a body-fitted computational mesh using 

r c 

= vortex core radius 

an uncoupled iterative solution. The interacting vortex 
wake is represented by an array of discrete vortices which. 

r 0 

= vortex core outer radius 

in turn, are represented by a finite-core model. The 
evolution of the interacting vortex wake is calculated by 

Re 

= Reynolds number, U«> 1/v 

Lagrangian techniques. The flow around circular and 
elliptic cylinders in the absence of an interacting vortex 

St 

= Strouhal number, nl/Uoo 

wake has been calculated. These results compare very 
well with other numerical results and with results obtained 

t 

= time 

from experiment and thereby demonstrate the accuracy of 
the viscous solution. The interaction of a simulated rotor 

U 

= free stream velocity 

wake with the flow about two-dimensional bodies, 
representing cross-sections of fuselage components, has 

u 

= velocity component in x -direction 

been calculated to address the vortex-interaction problem. 
The vortex interaction was calculated for the flow about a 

V 

= velocity component in y-direction 

circular and an elliptic cylinder at 45 and 90 degrees 
incidence. The results demonstrate the significant 

X 

= Cartesian coordinate 

variation in lift and drag on the two-dimensional bodies 
during the vortex interaction. 

y 

= Cartesian coordinate 

Symbols 

a 

= incidence angle 

Cd ■ drag coefficient 

P 

= transformation parameter 

Cf x = force coefficient in x -direction 

A 

- incremental value 

Cfy = force coefficient in y -direction 

r 

= circulation 

C/ = lift coefficient 


= coordinate in transformed space 

h£ = vector length in transformed space 

V 

= kinematic viscosity 

h|) = vector length in transformed space 

k 

= coordinate in transformed space 

J = transformation Jacobian 

1 = body length 

<D 

= vorticity 
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Subscripts 

vw = interacting vortex wake quantity 

0 = reference value 

1 = value at body surface 

oo = value at far Held 

Introduction 

The numerical prediction for the strong interaction 
between vortical wakes and the viscous flow field about 
bodies is of considerable importance in the design analysis 
of rotorcraft. The flow field surrounding a helicopter 
configuration is highly complex due to the lack of 
symmetry and the unsteady aspects of the flow. This 
unsteady, asymmetric flow field is complicated further in 
that the shed wake from the rotor blades interacts with 
other components of the aircraft (e.g., rotor blades, 
fuselage, tail boom, and tail rotor). In general, these are 
strong interactions in which the rotor wake flows onto or 
passes very close to the other components of the aircraft 
Vortex interactions occur at many flight attitudes (e.g., 
hover, descent, and low speed flight) and therefore have 
profound effects on the overall aerodynamic efficiency of 
the rotorcraft The effects of the vortex interactions are 
realized as increased vibratory loading on the fuselage, 
decreased payload capabilities, and increased noise. 

Because of the complexity of the problem and the limi- 
tations of the computer in modeling the entire flow field, 
the research effort addressing vortex interaction for 
rotorcraft has developed in two directions: blade-vortex 
interaction and rotor-fuselage interaction. Most of the 
vortex-interaction research for rotorcraft has been directed 
toward blade-vortex interaction in forward flight and, 
therefore, concentrated on the flow field near the rotor tip 
on the advancing side of the rotor. The flow field in this 
region is transonic and the rotor blade is at near zero 
incidence. Navier-Stokes calculations of blade-vortex 
interaction are being conducted with increasing efficiency 
(Refs 1, 2). These calculations model limited separation, 
due to the incidence of the blade and flow velocity. In 
addition, the blade-vortex interaction is generally not a 
direct interaction between the rotor blade and the rotor 
wake, but a very close passage. Aircraft descent is the 
exception when the rotor wake is convected through the 
rotor and directly interacts with the rotor blades. The flow 
beneath the rotor is quite different, with rotor-fuselage 
interactions present over a large portion of the fuselage 
length. 

The rotor-fuselage interaction has received far less 
attention than that for the blade-vortex interaction. These 
interactions, which exist for most flight configurations, are 


caused by the impingement of the convected rotor wake 
onto the fuselage components. An overview of the rotor- 
fuselage interactional aerodynamics can be found in 
Ref. 3. In contrast to the blade-vortex interaction, the flow 
field of the rotor-fuselage interaction is subsonic, highly 
separated, and vortex dominated. Portions of the fuselage 
are immersed in the rotor downwash, which is an unsteady 
flow with very large velocity gradients just inboard of the 
rotor tip. The fuselage is not streamlined in the direction 
of the downwash; as a consequence, separated wakes form 
at the surface owing to viscosity. These separated wakes 
are very unsteady due to the interaction of the impinging 
rotor wake. The result of the vortex interaction is the cou- 
pling of the vorticity produced at the surface of the fuse- 
lage with the interacting vortex wake. The flow becomes 
increasingly unsteady due to the mutual convection 
between the interacting vortex wake and the surface 
vorticity. The resulting flow field is very complex, and 
therefore difficult to model and predict. 

The analysis of rotor-fuselage interactions, as well as 
blade-vortex interactions, is complicated further by 
computer memory limitations. Currently available com- 
puter memory is inadequate for the three-dimensional 
meshes necessary to model the flow field surrounding a 
representative rotorcraft. Therefore, the computational 
mesh enclosing a body is generally concentrated near the 
body to resolve the flow at that location and considerably 
less dense away from the body. Owing to numerical 
diffusion, the coarse mesh away from the body is 
inadequate far accurate predictions for the evolution of the 
rotor wake surrounding a rotorcraft. In order to calculate 
the rotor-fuselage interaction, a method is needed to 
accurately predict both the rotor wake evolution away 
from the body on a coarse mesh and the separated viscous 
flow near the body surface. 

As a first step toward analyzing the full three-dimensional 
rotor-fuselage interaction, this paper presents a method to 
calculate, the unsteady, two-dimensional, incompressible 
interaction of an interacting vortex wake with the 
separated flow about bodies. The method solves the 
velocity/vorticity formulation of the Navier-Stokes 
equations. The computational mesh is concentrated near 
the body surface to resolve the boundary layer and is 
increasingly coarse further from the body. The interacting 
vortex write is modeled as an array of finite-core vortices 
(Refs 4, S). The core radius is variable and is independent 
of the mesh spacing. The finite-core model eliminates the 
numerical diffusion associated with the coarse mesh 
spacing away from the body and provides for the accurate 
prediction of the rotor wake evolution. The finite-core 
vortex model is accurate away from the body where the 
flow field is inviscid and dominated by the rotor wake. 
However, closer to the body, the flow is viscous and the 
convection of the rotor wake can no longer be considered 
inviscid. In this region, the interacting vortex wake 


2 


interacts with the viscous separated wake and can no 
longer be modeled by the finite-core vortex model. The 
vortex core is acted upon by shearing forces which will 
contort and possibly separate the vortex core into several 
separate vortices. As the vortex approaches the body, the 
finite-core vortex is distributed to the computational mesh 
and allowed to convect as part of the viscous solution. In 
this way the evolution of the rotor wake away from the 
body and the vortex interaction at the body are accurately 
simulated. In addition the present method for two- 
dimensional calculations is directly extendable to three- 
dimensional analysis. 

Results for the viscous flow about two-dimensional bodies 
have been calculated with and without an interacting 
vortex wake. The flow around circular and elliptic 
cylinders in the absence of an interacting vortex wake has 
been calculated in order to compare results of the current 
method with those of other numerical studies and 
experimental investigations. The intent of this paper is to 
provide comparisons between the flow about bodies with 
and without the presence of an interacting vortex wake. 
The simulated flow represents the down wash beneath a 
rotor. In this two-dimensional analysis, the body 
represents a cross-section of the tail boom or fuselage 
component, which may be circular or elliptical in 
geometry. 


Problem Formulation 

The vortex interaction problem is modeled by the 
velocity/vorticity form of the unsteady incompressible 
Navier-Stokes equations. The nondimensional governing 
equations in Cartesian coordinates are written: for the 
velocity. 


V^u = - 9to/9y , V^v = <ko/3x (la,b) 

and for the vorticity, 

coi + 9(uto)/9x + 9(vto)/9y= V 2 to/Re (2) 


where 

1 = the body length 

Re = Reynolds number based on body length 
Uoo = free stream velocity 
v = kinematic viscosity 

A body fitted computational mesh is used in the finite- 
difference solution to the governing equations. A 
hyperbolic mesh generator was used to generate the 
computational mesh, Ref 6. A portion of the computa- 
tional mesh is shown in Fig 1. The solution to the flow is 
obtained by solving the transformed Cartesian equations 
in the computational domain <£/n). 

The general transformation expression can be written, via 
the chain rule, as 

d( )/3£ = dx/9£(d( )/9x) +dy/3£(d( )/3y), and 
d( )fa\ = dx/dqO( )/dx) +dy/dr|(3( )/dy), or 


‘W 


9x/94 

9y/94 

9()/9x 

_9(V9n. 


9x/9q 

9y/9q_ 

_9(V9y 


The derivatives d( )/3x, and d( )/dy are solved for by 
inverting the coefficient matrix in the above equation. The 
Cartesian derivatives become 


3()/9x 

.dO0y. 


= W 


9y/9q 

-dx/dr\ 


-9y/94 

ax/94 


L9(yan 


where J is the Jacobian of the transformation. This relation 
is used to transform the derivatives of the governing 
equations to the generalized curvilinear coordinates. In the 
computational domain, the governing equations become 


where V 2 = 9^( )/fo 2 + 9^( )/9y2 


for the u-velocity component. 


where (x,y) are the Cartesian coordinates, R e is the 
Reynolds number, and t is the time. The variables (u,v) are 
the Cartesian components of the velocity and to, the 
vorticity, is defined by 

to = 9v/9x - 9u/9y (3) 

The nondimensional variables are written 


b n 2 (9 2 u/94 2 ) - 2p(9 2 u/949n) + h£ 2 (9 2 u/9n 2 ) 

+ Pdu/94 + Q9u/9q 

= -J[x£(9to/9n) - x^Oto/94)], (4a) 

for the v-velocity component. 


x = x’/l, y = y'/l, u = u'/Uoo, v = v'/Uoo 
to = to'/OWI), t = t7(l/Uoc), Re = Uool/v 
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hq^d^/t)^ 2 ) - 2P(d^v/3^dn) + h^d^/Sq 2 ) 
+ Pdv/3lj + Qdv/9q 


u = Uo«cos(a) 

(8a) 

v = Uoosin(a) 

(8b) 


= -J[yq(doV9£) - y£(9oVdn)] . (4b) 

and for the vorticity, 

©t + [yq(9uco/3£)- y£(du(D/<)q)]/J 
+ [x£(dvto/9q) - xq(9v(o/9£)]/J 

= [hn 2 ^ 2 ^ 2 ) - 2p(d 2 co/3^dn) + h^2(d 2 avan 2 ) 

+ P3aV9£ + Q3oV9T|l/(ReJ 2 ) (5) 

where (£,q) are the transformed coordinates, h£ and bq 
are the vector lengths, 

h£ = (xjj 2 + y£ 2 ) 1/2 and 
bq = (*q 2 + yq 2 ) 1/2 

with 

P = x^xq + yijyq, 

J = x^yq - xqy^, 

P = [y^hq 2 *^ - 2Px^q + hq^qq) 

- x^fbq^^ - 2Py^q + hq^qq)]/!, 

Q = [-yriOtq^ - 2pxjjq + hq^qq) 


and 

(O = 0 (9) 

in the potential flow region upstream of the body, i.e., for 
q 1 < T| < t \2 of Fig 1, with q increasing in the counter- 
clockwise direction. 

In the region downstream of the body, the velocity at the 
far field boundary is calculated from Eqs 4a and 4b and 
the vorticity is calculated from Eq 5 by considering the 
inertia terms only and neglecting the viscous terms. 
Therefore, the vorticity at the far field boundary, for 
ti? < n < r imyr and qmin < q <qi is calculated from 

cot. + [yq(duto/34) - y£(duoVdq)]/J 

+ [x^(dvoV9q) - xq(3v&V9^)]/I = 0. (10) 

The surface pressure is obtained by integrating the tan- 
gential component of the Navier-Stokes equation along 
the surface of the body. At the surface 

pq = -[(po>q - hq20GV^)]/(ReJ) (11) 

or 

PI =Po- 1/ReJ [p®q - hq 2 (dciV9£)/Hdq (12) 

% 


+ xqfhq^^ - 2Py£q + hq^qqXl/I, 
and 

x£ = 9x/9£, xq = 3xy^n, y£ = dy/3£, yq = 9y/9q, 
x££ = 9 2 x/9£ 2 * xqq = 9 2 x/9q 2 * y££ = 9 2 y/9£ 2 * 
yqq = dhftn 2 

The boundary conditions for Eqs 4a, 4b, and S at the body 
surface are calculated from the no-slip condition as 

u = 0, v = 0 (6) 

and the definition of vorticity for to, 

to = [yq9v/3£ + xq9u/9£]/J (7) 


where pi is the pressure on the surface and po is the 
reference pressure on the surface. 

The reference pressure is obtained by the line integral of 
the momentum equation from the far field boundary to the 
body surface along a constant q line. The reference 
pressure is calculated from 

r ^oo 

Po= I (9(x^u + y£v)/9t - (x£vo) - y^uco) 

\ 

+ h£ 2 0(fl0n) - P0oV3S)/(ReJM + 1/2 (13) 

The force components acting on the body are calculated 
from the integration of the pressure and shearing forces at 
the body surface. The force coefficients are 


The boundary conditions at the far field boundary are 
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Numerical Method 


f ^max 

Cfe=-J (piyil)<hl 

^min 

f ^max 

- (1/Re) J (wixr,^ (14) 


C ^max 

Cfy=-J (Pix^dii 

^min 

C ^max 

-(1/Re) J (wiyTi)dn (15) 

^min 

where CD] is the vorticity at the body. 


The lift and drag coefficients are 


Q = - Cf x [sin(a)) + Cf y (cos(a)] 

(16) 

Cd = Cf x [cos(a)) + Cfy(sin(a)] 

(17) 

where a is the incidence of the body. 


The interacting vortex wake is modeled as an array of 
finite-core vortices. A third-order-polynomial distribution 
of vorticity about each vortex position is used to represent 
the vortex core, Ref 4. This distribution of vorticity is 

tOvwW = f /(2nr 0 ^)[ 20/3(r3/r Q 3) 


- IOO^/tq 2 ) + 10/3] 

(18) 


where a> V w represents the vorticity of the interacting 
vortex wake and r 0 is the outer radius of the vortex dis- 
tribution. The core radius, r c , the radial position of 
maximum tangential velocity, is equal to 0.6376 r 0 for this 
finite -core vortex model. The influence of the interacting 
vortical wake is included in Eqs (4a) and (4b) through the 
4 and i) components of -( 7 x ffiyw)- The interacting 
vortex wake is modeled as an array of finite-core vortices 
prior to the interaction. However, when the interacting 
vortex wake is less than a prescribed distance from the 
body surface, the finite-core model is no longer adequate 
to model the wake. At this point, the interacting vortex 
wake is distributed to the computational mesh and 
included as part of the viscous solution. In this way the 
vortex interaction is modeled by the governing equations. 


The governing equations are solved by standard finite 
difference methods on the body-fitted computational 
mesh. The mesh is concentrated near the body to resolve 
the flow at that location. The solution is started 
impulsively from rest The impulsive start is represented 
by a potential velocity distribution in the flow field with 
zero velocity and a sheet of vorticity at the body surface. 

The vorticity transport equation, Eq (S), is solved using an 
alternating direction implicit method (ADI) (Ref 7). The 
ADI method solves the time-dependent vorticity transport 
equation by separating the equation into two implicit 
equations representing two fractional steps for a single 
time step. This method results in the solution of scalar 
tridiagonal matrices for each fractional step. 

The vorticity at the body surface in Eq (S) is calculated 
from Eq (7), The velocity, Eqs (4a) and (4b), are cast in a 
time-dependent form and solved by the ADI method. The 
influence of the interacting vortex wake is included in the 
velocity equations through the negative curl of the 
vorticity distribution in the vortex core, 
-<V x (Ovw)- The vorticity distribution in the vortex core is 
defined in Eq (18). Outside the outer radius the vorticity is 
zero. 

The solution to the time-dependent problem is calculated 
by solving the vorticity transport equation and the velocity 
equations in an uncoupled manner. In this way, the 
vorticity transport equation is solved first, after which 
each velocity equation is solved in turn. This solution 
procedure continues until a converged solution is obtained 
for the current time step. The influence of the interacting 
vortex wake is considered constant at each time step 
during this iteration process. At convergence, the vorticity, 
flow field velocity, and interacting vortex wake position 
are updated and the solution advances in time. 

Computed Results 

The method described above has been applied to the flow 
about a circular cylinder and a 25% elliptic cylinder for a 
Reynolds number of 3000. The computational grid is 
dimensioned 81 x 241 with the far field radius roo = 20 for 
all calculations. The time is nondimensionalized with 
respect to the length of the body, in this case the diameter 
of the cylinder. Therefore, a time interval of At = 1 .0, 
corresponds to the time required to travel a distance equal 
to the cylinder diameter at a velocity equal to that of the 
free stream velocity. The flow is calculated with and 
without the presence of the interacting vortex wake. 

The flows about a circular cylinder and an elliptic cylinder 
at 45 and 90 degrees incidence are calculated without the 
interacting vortex wake and are included for comparison 
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with the results of the vortex interaction. The incidence 
angle is measured with respect to the vertical axis as 
shown in Fig 1 . Comparisons of the results for the circular 
cylinder without the interacting vortex wake can be found 
in Ref 8. In this reference comparisons are made between 
the results of the present method, the experimental results 
of Ref 9 and the numerical results of Ref 10. In Ref 8, 
comparisons are made between the current method and the 
results of Ref 9 for the surface vorticity and were shown 
to be excellent In addition, the velocity profile in the 
wake of the circular cylinder for the current method was 
compared with the results of Refs 8 and 9. The 
comparisons were shown to be excellent 

The flow around a circular cylinder and a 25% elliptic 
cylinder is calculated during the vortex interaction. The 
flow about the elliptic cylinder during the vortex inter- 
action was calculated for 45 and 90 degrees incidence 
with respect to the vertical axis. The flow about the 
elliptical cylinder at 45 degrees incidence was calculated 
for a vortex strength of each vortex in the interacting 
vortex wake of r = —1.0 and r = 1.0. The flow about the 
circular cylinder and the elliptic cylinder at 90 degrees 
incidence is only calculated for a vortex strength of each 
vortex in the interacting vortex wake of T = 1.0 due to 
symmetry. A vortex of strength T = -1 represents a 
clockwise rotation, while a vortex of strength r = 1 
represents a counterclockwise rotation. These vortex 
interactions were considered to evaluate general vortex 
interactions and because no specific experimental data 
representative of this flow exists. 

No Vortex Interaction Case 

Circular Cylinder. The variation of the lift and drag 
coefficients with time for the flow about the circular 
cylinder are shown in Fig 2. The effect of the secondary 
vortices on the drag coefficient can be seen as the small 
variations in the curve. The development of the separated 
vortices remain symmetric up to t = 40.0 and, therefore, 
the lift coefficient is zero. Comparisons between these 
results and the results of other investigations can be found 
in Ref 8. The development of the streamlines of this flow 
are shown in Fig 3. The figure demonstrates the 
emergence of secondary vortices downstream of the 
cylinder. These results will be compared with the results 
during a vortex interaction presented later in this paper. 

Elliptic Cylinder. The time evolution of the lift and drag 
coefficients far the flow about a 25% elliptic cylinder are 
shown in Fig 4 far 45 degrees incidence. For t greater than 
t = 20.0, the Strouhal number (St = nl/Uoo) for this flow is 
equal to 0.174 based on the projected body length. The 
development in time of the streamlines around the elliptic 
cylinder are shown in Fig 5 for 45 degrees incidence. 


Variations in the lift and drag coefficients with time for 
the flow around an elliptic cylinder at 90 degrees inci- 
dence are shown in Fig 6. For this flow the separated 
vortices do not remain symmetric and begin to oscillate at 
t * 5.0. This is indicated in the nonzero lift coefficient and 
is shown in the development of the streamlines in Fig 7. 
The Strouhal number of this flow is equal to 0.214. 

Vortex Interaction Case 

The interaction of a vortex wake with the flow about a 
circular cylinder and the flow about a 25% elliptic 
cylinder are calculated. The interaction of the vortex wake 
with the flow about the elliptic cylinder was calculated for 
45 and 90 degrees incidence with respect to the free 
stream. 

Each interaction modeled the interacting vortex wake as 
an array of vortices inserted into the flow field at constant 
time intervals. This represents the periodic passage of a 
rotor blade and the associated tip vortex, thus generating 
the rotor’s vortical wake structure. The orientation of the 
interacting vortex wake with respect to the elliptic 
cylinder is depicted in Fig 8. The free stream velocity is 
directed along the negative y-axis and represents the 
downwash due to a rotor. The vortex interaction was 
started at t = 0.25, which was long enough after the 
impulsive start as not to adversely affect the initial 
development of the flow. At the start of the interaction, 
five vortices are inserted in the flow field, as shown in Fig 
8. Then, as the interaction proceeds, an additional vortex 
is repeatedly inserted in the flow field at the uppermost 
position, y = 6.5, at a constant time interval. At = 1.0. 
Therefore, die interacting vortex wake will represent the 
periodic rotor wake during the entire vortex interaction. 

Circular Cylinder , r = 1.0. The variation in the lift 
coefficient and drag coefficient during the vortex inter- 
action are shown in Figs 9 and 10, respectively. The 
immediate effect of the interaction is to rapidly increase 
the lift on the cylinder. Fig 9a, and slightly reduce the 
drag. Fig 10a. This increased lift is the result of the 
inrhirf/t angle of attack of the cylinder due to the presence 
of the interacting vortex wake. The lift on the cylinder in 
the absence of the interacting vortex wake was zero. As 
the interaction progresses, the effect of the periodic 
passage of the interacting vortex wake can be seen in the 
va riatio n in the lift on the cylinder. The lift coefficient 
rises rapidly and then is reduced quickly in a sawtooth 
pattern representative of a vortex interaction. The 
separation in these lift peaks is approximately t = 1.0. This 
corresponds to the natural convection of the vortices at the 
free stream velocity, U<x> = 1.0. The effect of the vortex 
interaction on the drag. Fig 10, is less pronounced than 
that for the lift. Fig 9. At t = 5.0 a strong vortex interaction 
begins to occur. This strong vortex interaction is the result 
of the initial vortex trajectory taken by the uppermost 
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vortex at the start of the interaction. As the interaction 
starts, the uppermost vortex trajectory is acted upon by the 
free stream velocity and the influence of the vortices 
beneath this position. Because there are no vortices above 
this position, the vortex is convected in a direction down 
and to the left. This is indicated in the vortex trajectories 
of Fig 11 and the development of the streamlines during 
the interaction, as shown in Fig 12. In Fig 12e at t = S.0, a 
pair of vortices can be seen just above the cylinder 
surface. The vortex to the left is the vortex that started at 
the uppermost position and is being convected into the 
cylinder due to the effect of the vortex to the right. After 
the strong interaction, the interacting vortex wake has 
regained the the constant vortex separation similar to the 
initial distribution of the interacting vortices. However, 
the vortices of the interacting vortex wake are now 
convected along a curved trajectory to the left of the initial 
distribution, as shown in Fig 1 lb. The lift returns to the 
sawtooth distribution prior to the strong interaction. Fig 9, 
t = 6.0. For time t > 20.0, a pattern of periodic vortex 
shedding from the cylinder is shown in the lift. Fig 9b, and 
the drag. Fig 10b. 

The vortex trajectories during the interaction are shown in 
Fig 11. The trajectories of the five vortices initially 
inserted into the flow at t = 0.2S are shown in Fig 11a and 
the trajectories of the vortices inserted in the flow at y = 
6.S for the duration of the interaction are shown in Fig 
1 lb. The figure indicates that the initial vortex that starts 
at the uppermost position, y = 6.5, is convected in a 
direction down and to the left. It is this vortex that 
generates the strong vortex interaction because it is 
convected nearest to the cylinder surface owing to the 
effects of the other vortices. Once the initial vortex at the 
uppermost position has passed the cylinder, the interacting 
vortices of the vortex wake convect in a regular pattern 
along the curved trajectory of Fig lib. 

Development of the streamlines for the vortex interaction 
are shown in Fig 12. Initially, Fig 12a, t * 1.0, the vortices 
are aligned in a pattern similar to that of the initial vortex 
distribution. Later, at t = 2.0, Fig 12b, a vortex has been 
distributed to the mesh and is shown as the local curvature 
in the streamlines just above the cylinder surface. At t = 
4.0, the vortex that started at the uppermost position, y « 
6.5, can be seen near the stagnation streamline and is 
being convected at an increased rate as to overtake the 
adjacent vortex. This increased descent is due to the 
influence of the adjacent vortex. The vortex which started 
at the uppermost position passes closest to the cylinder 
surface at t = 6.0 and is the closest approach of any vortex 
during the vortex interaction. Later on, the vortices of the 
interacting vortex wake are very regularly spaced and are 
being convected along the curved trajectory toward the 
cylinder surface. 


Elliptic Cylinder. Several vortex interactions are con- 
sidered for flow about a 25% elliptic cylinder The dif- 
ferences in the vortex interactions are the result of varying 
the incidence angle and the rotation sense of the vortices 
that comprise the interacting vortex wake. The vortex 
interaction is considered for 45 and 90 degrees incidence. 
Due to symmetry, a single interaction is considered for the 
flow at 90 degrees incidence. At 45 degrees incidence, 
two interactions are considered, the difference being the 
rotation sense of the vortices in the interacting vortex 
wake. The three interactions are described below. 

First Interaction, T = -1.0, a = 45°: The vortex interaction 
is started att = 0.25, as before. The variation in the lift and 
drag coefficients during the vortex interaction are shown 
in Figs 13 and 14, respectively. The initial effect of the 
vortex interaction is to rapidly reduce the lift on the 
elliptic cylinder. The is due to the movement of the 
stagnation streamline and therefore, a reduction in the 
induced incidence of the elliptic cylinder. A corresponding 
reduction in the drag is seen in Fig 14. The reductions in 
the lift and drag are quickly reversed as the vortex 
interaction progresses, t = 2.0, and the sawtooth signature 
in the lift and drag develop. The sawtooth pattern is 
representative of the passage of the periodic interacting 
vortex wake over the elliptic cylinder. As in the previous 
vortex interaction for the circular cylinder, the vortex that 
starts at the uppermost position, y = 6.5, convects in a 
direction so as to pass closest to the surface of the elliptic 
cylinder. During this interaction, the uppermost vortex is 
convected in a direction down and to the right, which is 
again toward the stagnation streamline and is caused by 
the effects of the interacting vortices below the initial 
vortex position. This is clearly seen in the development of 
the vortex trajectories of Fig 15 and the development of 
the streamlines shown in Fig 16. The strong vortex 
interaction seems to occur at t = 6. At this time a change 
in the lift. Fig 13, and the drag. Fig 14, is observed. A 
spike in the loading before t = 4.5 and after t = 8.5 is not 
present in the lift and drag distribution at t = 6.0 due to the 
strong vortex interaction. At t = 8.0, the lift. Fig 13, and 
the drag. Fig 14, return to the sawtooth distribution prior 
to the strong interaction. As in the previous interaction, 
the interacting vortex wake has regained the constant 
vortex separation, as shown in Fig 16, similar to the initial 
distribution of the interacting vortices, and the vortices are 
being convected along a curved trajectory. The curved 
trajectory is now to the right of the initial distribution, as 
shown in Fig 15b. For t > 20.0 the loading on the elliptic 
cylinder is periodic. The effect of the vortex interaction 
can be seen as the “spikes** in the loading. 

The vortex trajectories during the interaction are shown in 
Fig 15, and indicate that the initial vortex which starts at 
the uppermost position, y = 6.5 is convected in a direction 
down and to the right, as shown in Fig 15a. It is this 
vortex that generates the strong vortex interaction because 
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it is converted along trajectories close to the stagnation 
streamline of the cylinder. This trajectory provides the 
closest approach to the surface of the elliptic cylinder as 
shown in the figure. After this interaction, the vortices are 
converted along the curved trajectory shown in Fig 15b. 

Development of the streamlines for the vortex interaction 
are shown in Fig 16. At the early developments of the 
interaction. Figs 16a, 16b, and 16c (t = 1.0, t = 2.0, and t * 
3.0), the vortices are aligned in a pattern similar to that of 
the initial vortex distribution. When t = 4.0 (Fig 16d) the 
vortex started at the uppermost position can be seen 
closest to the stagnation streamline. It is this vortex that 
will make the closest approach to the cylinder and 
generate the strong vortex interaction. At t = 5.0, the effect 
of the adjacent vortex to the left of the vortex started at the 
uppermost position is shown to convert this vortex toward 
the cylinder. The interacting vortex wake has a profound 
effect on the development of the separated wake below the 
elliptic cylinder. This is most significant at t = 7.0, 8.0, 
and 9.0, where three separated vortices are in close 
proximity to the lower surface of the elliptic cylinder. 
During the developing flow in the absence of the 
interacting vortex wake (Fig 5), the shed vortices oscillate 
from leading to trailing edge, thereby generating a vortex 
street in the wake of the cylinder. This vortex street is 
altoed significantly during the vortex interaction. 

Second Interaction, T = 1.0, a = 45°: This vortex inter- 
action is calculated because the flow about the elliptic 
cylinder at 45 degrees incidence is not symmetric with 
respect to the rotation direction of the vortices in the 
interacting vortex wake. The vortex interaction is again 
started at t » 0.25. The vortex strength is now equal to T = 
1.0, for each vortex in the interacting vortex wake. This 
corresponds to a counter-clockwise rotating vortex of 
magnitude r = 1.0. The variation in the lift and drag 
coefficients during the vortex interaction are shown in 
Figs 17 and 18, respectively. During this vortex interac- 
tion, the effect of the interacting vortex wake is to sig- 
nificantly increase the lift and drag on the cylinder. This is 
due to the rotation of the stagnation streamline which 
increases the induced incidence of the cylinder. At 
approximately t = 1.5, die sawtooth pattern in the lift and 
drag can be seen. The pattern is less pronounced during 
this interaction than the previous interaction because the 
effect of the interacting vortex wake is less significant in 
comparison with the vorticity generated at the surface of 
the elliptic cylinder. The strength of the shed vorticity at 
the leading edge of the cylinder increases with incidence 
and will dominate the flow near this location. The 
sawtooth pattern is fairly consistent for the lift. Fig 17, and 
for the drag, Fig 18, up to a time oft= 5.5. At this timet = 
5.5, the strong vertex interaction is occurring, although 
not clear from the lift and drag. This strong interaction is 
again the result of the convection of the uppermost vertex 
due to the influence of the other interacting vortices and 


the closest approach of this vertex to the surface of the 
cylinder. This can be seen in the vortex trajectories of 
Fig 19 and the streamline contours of Fig 20. The effect of 
the interacting vortices during the interaction is to first 
accelerate the flow around the leading edge of the elliptic 
cylinder and to then retard this flow. Figs 20e and 20f. 
This has a significant effect on the separated vortex wake 
generated at that location and, therefore, the vortex wake 
below the cylinder. 

The variation in the drag coefficient on the cylinder during 
the interaction is shown in Fig 18. As in the previous 
interaction. Figs 13 and 14, and unlike the interaction of 
the vortex wake with the circular cylinder, Figs 9 and 10, 
the lift. Fig 17, and the drag, Fig 18, show similar trends 
during this interaction. At the early stages of the 
interaction, t = 1.0, the drag is sharply increased due to the 
interaction, after which the sawtooth pattern is present, t = 
1.5 tot = 5.0. At t = 5.0, the drag is significantly increased 
due to the strong vertex interaction as was seen in the lift. 
Fig 17. Lam, the sawtooth pattern in the drag reemerges, t 
> 8.0, Fig 18. Later, Fig 18b, the loading on the cylinder 
becomes periodic with a Strouhal number similar to that 
for the flow with no vertex interaction. 

The vertex trajectories during the interaction are shown in 
Fig 19 and indicate that the vortex which starts at the 
uppermost position, y = 6.5, is converted in a direction 
down and to the left due to the influence of the other 
interacting vertices. It is this initial vortex that generates 
the strong vertex interaction because it is converted along 
trajectories close to the stagnation streamline of the 
cylinder and passes closest to the cylinder surface. After 
this vertex passes the cylinder, the interacting vortices are 
converted along a curved trajectory. Fig 19b, to the left of 
the initial distribution. 

Development of the streamlines for the vertex interaction 
are shown in Fig 20. At the early developments of the 
interaction. Figs 20a, 20b, and 20c (t = 1.0, t = 2.0, and t = 
3.0), the vertices are aligned in a pattern similar to that of 
the initial vortex distribution and have significantly 
increased the induced incidence of the elliptic cylinder. 
When t = 4.0 (Fig 20d) the vortex that started at the 
uppermost petition can be seen as the vortex to the left in 
the figure. Att= 5.0 and t = 6.0, this vortex is close to the 
surface of the elliptic cylinder and will pass closest to the 
cylinder. It is this vortex that contributes to the strong 
vortex interaction. 

Third Interaction, T = 1.0, a = 90°: The variation in the 
lift and drag coefficients during the vortex interaction are 
shown in Figs 21 and 22, respectively. The effect of the 
interaction is to first increase the lift on the elliptic 
cylinder, up to t = 6.0, at which time a strong vortex 
interaction occurs. The sawtooth pattern in the lift, 
indicative of vortex passage, is again present. In 
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comparing the lift during this interaction. Fig 21, with the 
lift during the vortex interaction with the circular cylinder. 
Fig 9, the current interaction generates a less pronounced 
effect on the lift This is due to the curvature of the 
respective cylinders and the effect of the interacting vortex 
wake on these surfaces. The variation in the drag during 
this interaction is shown in Fig 22. The drag is initially 
increased, but only slightly when compared with the drag 
for no interaction as shown on the figure. Later, large 
spikes in the drag distribution are present at t = 2.5, 3.5, 
4.5, and 5.5, representing the passage of the interacting 
vortices. These spikes in the drag distribution are due to 
the influence of each vortex in the interacting vortex 
wake. As each vortex approaches the elliptic cylinder, the 
flow at the upper surface is accelerated resulting in a 
corresponding reduction in drag. The drag during the 
interaction nearly returns to the value for no interaction. 
The strong vortex interaction at t = 6.5 disrupts this effect 
and suppresses the spike in the drag. At t = 8.0, the strong 
vortex interaction is over and the drag. Fig 22, begins to 
return to the distribution shown prior to the interaction. 
For t > 20.0, the flow becomes periodic with a Strouhal 
number similar to that for the no-interaction flow. 

The vortex trajectories during the interaction are shown in 
Fig 23 and indicates that the initial vortex which starts at 
the uppermost position, y = 6.5, is convected in a direction 
down and to the left due to the influence of the other 
interacting vortices. This vortex passes close to the upper 
surface of the cylinder, disrupting the flow at that location. 
Fig 23a. The interacting vortices are eventually convected 
along a curved trajectory. Fig 23b, to the left of the initial 
vortex distribution. Closer to the cylinder, the vortices do 
not follow a single trajectory, similar to that for the 
interaction with the flow about a circular cylinder. Fig 11. 
Instead, the vortex trajectories are fanned out near the 
surface. 

Development of die streamlines for the vortex interaction 
are shown in Fig 24. At the early developments of the 
interaction. Figs 24a, 24b, and 24c (t= 1.0, t = 2.0, and t = 
3.0), the vortices are again aligned in a pattern similar to 
that of the initial vortex distribution. The effect of the 
interacting vortices is to increase the upper surface 
circulation around the right side of the cylinder. This 
produces the increase in lift. Fig 21, and the corre- 
sponding spikes in the drag. Fig 22 for the passage of each 
interacting vortex. The strong vortex interaction is shown 
for times t = 6.0, Fig 24f when the initial vortex at the 
uppermost position passes closest to the cylinder. Figure 
24f clearly shows the effect of the interacting vortex wake 
on the shed wake below the cylinder when these results 
are compared with the flow for no interaction. Fig 7. 


Conclusions 

A method has been developed to model the two- 
dimensional interaction between an interacting vortex 
wake represented by a finite-core model and the viscous 
flow around arbitrary bodies. The method solves for the 
flow field velocities on a body-fitted computational mesh 
using finite-difference techniques. The viscous flow field 
of the two-dimensional body is calculated on an Eulerian 
grid via the velocity/vorticity formulation of the Navier- 
Stokes equations. 

The method demonstrates that the finite-core vortex model 
is accurate in convecting the interacting vortex wake away 
from the body. This will significantly reduce the need for 
a fine computational mesh away from the body to resolve 
the interacting vortex wake. Near the body, the interacting 
vortex wake is distributed to the computational mesh to 
provide accurate vortex interaction with the flow field. 

A simulation of a rotor wake interaction with the flow 
about a circular cylinder and a 25% elliptic cylinder at 45 
and 90 degrees incidence was shown for Reynolds number 
3000. The simulation was considered to address the 
interaction of a rotor wake with the cross-section of the 
tail boom and empennage components. Results of the 
vortex interaction indicate that the presence of a 
interacting vortex wake has a profound effect on the flow 
field and therefore, the loading on the circular and elliptic 
cylinders. The calculation of vortex interactions with the 
flow about a circular and an elliptic cylinder were 
considered to evaluate general vortex interactions. 

The lift and drag on the circular cylinder and the elliptic 
cylinder at 90 degrees incidence were significantly 
increased due to the interacting vortex wake. The lift on 
the circular cylinder during the interaction was increased 
significantly and the drag was slightly reduced from that 
for no interaction due to the effect of the interacting vortex 
wake on the circulation around the cylinder. The drag on 
the elliptic cylinder at 90 degrees incidence was sig- 
nificantly increased due to the vortex interaction. Spikes 
in the drag distribution, representing the passage of the 
vortices in the interacting vortex wake were shown. The 
magnitude of these spikes are significant with respect to 
the drag rise due to the vortex interaction. 

Large variations in the loading of the elliptic cylinder at 
45 degrees incidence were shown for changes in the 
rotation direction of the vortices in the interacting vortex 
wake. When the vortices were rotating in a clockwise 
direction, the lift and drag on the elliptic cylinder were 
significantly reduced during the vortex interaction. 
Conversely, when the vortices were rotating in a counter- 
clockwise direction, the lift and drag on the elliptic 
cylinder were significantly increased during the vortex 
interaction. These variations are important in 


9 



und erstanding the interaction of the rotor wake with the 
flow about empennage components. 

The present two-dimensional method has demonstrated 
the ability to efficiently compute the two-dimensional 
interaction of finite-core vortex wakes with the flow about 
arbitrary bodies. The method is currently being extended 
to three-dimensional analysis. 
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Fig 2 Lift and drag coefficients for circular cylinder. 
Re ■ 3000, a ■ 0°. 
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Fig 3 Streamline contours for circular cylinder, R e = 3000, a = 0°. (a) t = 1.0, (b) t = 2.0, (c) t = 
(d) t = 4.0, (e) t = 5.0, (0 1 = 6.0. 
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Fig 3 Concluded, (g) t = 7.0, (h) t = 8.0, (i) t = 9.0, (j) t = 10.0. 



Fig 4 Lift and drag coefficients for an elliptic cylinder, R e = 3000, a = 45°. 
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Fig 5 Streamline contours for an elliptic cylinder, R e = 3000, a = 45°. (a) t = 1.0, (b) t = 2.0, (c) t « 3.0, 
(d) t = 4.0, (e) t = 5.0, (0 1 = 6.0. 
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Fig 7 Streamline contours for an elliptic cylinder, Re - 3000, a = 90°. (a) t - 1.0, (b) t - 2.0, (c) t = 3.0, 
(d) 1 = 4.0, (e) t= 5.0, (0 1 - 6.0. 
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Fig 8 - Schematic of vortex interaction. 
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Fig 11 Vortex trajectories for vortex interaction of vortex strength T = 1.0, r<) = 0.1, circular cylinder. Re = 3000, 
a = 0°. (a) initial vortices, t = 0.25, (b) additional vortices. At = 1.0. 
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Fig 12 Streamline contours for vortex interaction of vortex strength T = 1.0, r D = 0.1, circular cylinder. Re = 3000, 
o = 0°. (a) t = 1.0, (b) t = 2.0, (c) t = 3.0, (d) l = 4.0, (e) t = 5.0, (f) t = 6.0. 
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Fig 12 Concluded, (g) t = 7.0, (h) t = 8.0, (i) t = 9.0, (j) t = 10.0. 
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Fig 13 Lift coefficient for vortex interaction of vortex 
strength T = -1.0, r 0 = 0.1, elliptic cylinder. Re = 3000, 
a = 45°. (a) t£ 10.0, (b)t< 40.0. 


Fig 14 Drag coefficient for vortex interaction of 
vortex strength T » -1.0, r 0 = 0.1, elliptic cylinder, 
Re = 3000, a = 45°. (a) t £ 10.0, (b) t < 40.0. 
















Fig 15 Vortex trajectories for vortex interaction of vortex strength T - -1.0, r 0 = 0.1, elliptic cylinder. Re = 3000, 
a = 45°. (a) initial vortices, t = 0.25, (b) additional vortices, At = 1.0. 
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Fig 16 Streamline contours for vortex interaction of vortex strength r = -1.0, r 0 = 0.1, elliptic cylinder. Re = 3000, 
a = 45°. (a) t = 1.0, (b) t = 2.0. (c) t = 3.0, (d) t = 4.0, (e) t = 5.0, (f) t = 6.0. 
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Fig 16 Concluded, (g) t = 7.0, (h) t = 8.0, (i) t = 9.0, <j) t = 10.0. 
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Fig 17 Lift coefficient for vortex interaction of vortex 
strength T = 1.0, r 0 = 0.1, elliptic cylinder, Re = 3000, 
a = 45°. (a) t < 10.0, (b) t < 40.0. 


Fig 18 Drag coefficient for vortex interaction of 
vortex strength T = 1.0, r D = 0.1, elliptic cylinder. Re = 
3000, a = 45°. (a) t £ 10.0, (b) t < 40.0. 














Fig 19 Vortex trajectories for vortex interaction of vortex strength r = 1.0, r 0 = 0.1, elliptic cylinder. Re - 3000, 
a = 45°. (a) initial vortices, t = 0.25 (b) additional vortices, At = 1.0. 
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Fig 20 Streamline contours for vortex interaction of vortex strength I~ = 1.0, r G = 0.1, elliptic cylinder. Re = 3000, 
a = 45°. (a) l = 1 .0, (b) t = 2.0, (c) t = 3.0, (d) t = 4.0, (e) t = 5.0, (0 1 = 6.0. 
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Fig 20 Concluded, (g) t = 7.0, (h) t = 8.0, (i) t = 9.0, (j) t = 10.0. 
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Fig 21 Lift coefficient for vortex interaction of vortex Fig 22 Drag coefficient for vortex interaction of 

strength T = 1.0, r 0 = 0.1, elliptic cylinder. Re = 3000, vortex strength T = 1.0, r 0 = 0.1. elliptic cylinder, Re = 

a = 90°. (a) t < 10.0, (b) t < 40.0. 3000, a = 90°. (a) t < 10.0, (b) t < 40.0. 
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Fig 23 Vortex trajectories for vortex interaction of vortex strength r= 1.0, r 0 = 0.1, elliptic cylinder, Re = 3000, 
a = 90°. (a) initial vortices, t = 0.25, (b) additional vortices. At = 1.0. 
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Fig 24 Streamline contours for vortex interaction of vortex strength T = 1.0, r 0 = 0.1, elliptic cylinder, R e = 3000, 
a = 90°. (a) t = 1.0, (b) t = 2.0, (c) t = 3.0, (d) t = 4.0, (e) t = 5.0, (0 l = 6.0. 
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Fig 24 Concluded, (g) t = 7.0, (h) t = 8.0, (i) t = 9.0, (j) t = 10.0. 
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